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Abstract. - In biomimetic and biological systems, interactions between surfaces are often medi- 
ated by adhesive molecules, nanoparticles, or colloids dispersed in the surrounding solution. We 
present here a general, statistical-mechanical model for two surfaces that interact via adhesive 
particles. The effective, particle-mediated interaction potential of the surfaces is obtained by inte- 
grating over the particles' degrees of freedom in the partition function. Interestingly, the effective 
adhesion energy of the surfaces exhibits a maximum at intermediate particle concentrations, and 
is considerably smaller both at low and high concentrations. The effective adhesion energy corre- 
sponds to a minimum in the interaction potential at surface separations slightly larger than the 
particle diameter, while a secondary minimum at surface contact reflects depletion interactions. 
Our results can be generalized to surfaces with specific receptors for solute particles, and have 
direct implications for the adhesion of biomembranes and for phase transitions in colloidal systems. 
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Introduction. — The adjustment of surface interac- 
tions is crucial for controlling the phase behavior of col- 
loidal systems [1] and the adhesiveness of biological cells 
and membranes [2]. These interactions are often domi- 
nated by the composition of the surfaces, which may be 
charged, hydrophobic, etc. In some systems, the interac- 
tions are also strongly affected by molecules or particles 
in the surrounding medium. The concentration of these 
particles is an additional control parameter for the surface 
interactions, a parameter that is often easier to adjust than 
the surface composition, and can be varied over a wider 
range than external parameters such as temperature. 

On the one hand, non-adhesive particles can induce at- 
tractive 'depletion' interactions between surfaces, because 
close contact of the surfaces reduces the excluded volume 
for the particles [3,4]. On the other hand, adhesive parti- 
cles can directly bind two surfaces together. For example, 
the surface interactions and colloidal phase of membrane- 
coated silica beads [5, 6] and the force between membrane- 
coated mica surfaces [7] have been altered by adding solu- 
ble, adhesive proteins. Multivalent ions such as chromium 
can induce the adhesion of lipid membranes, presumably 
by crosslinking the polar headgroups of lipids in apposing 
membranes [8]. Linker proteins that interconnect mem- 
brane receptors are known to assist biomembrane adhe- 
sion [2] , which has also been utilized in biomimetic exper- 
iments [9, 10]. 




Fig. 1: Two surfaces in contact with a solution of adhesive 
particles. A particle can bind the two surfaces together for 
surface separations slightly larger than the particle diameter 
(particle on the right). At large separations, the particles can 
only bind to one of the surfaces (particles on the left). 



Adhesive particles can bind two surfaces together if the 
separation of the surfaces is equal to or slightly larger than 
the particle diameter, see fig. 1. At larger separations, the 
particles can only bind to one of the surfaces. In this letter, 
we consider particles that exhibit short-ranged, attractive 
interactions with the surfaces, and repulsive hard-sphere 
interactions with each other. Our central result is that the 
effective, particle-mediated adhesion energy of the surfaces 
is given by 
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with a dimensionless coefficient a for small bulk volume 



p-1 



Bartosz Rozycki, Reinhard Lipowsky and Thomas R. Weikl 




OOOOOOCX) 



(XXXXXXX) 



CXXXXXXX3 



binding separation 



Fig. 2: Effective adhesion energy Ud of the surfaces, given in eq. (1), as a function of the particle bulk volume fraction 
for the binding energy U = 8T and q — 0.25. The effective adhesion energy is maximal at the optimal bulk volume fraction 
(f>* ~ e'^^"^ /q ~ 1.34 • 10~^. At the optimal volume fraction, the particle coverage of two planar parallel surfaces is 50% for 
large separations, and almost 100% for small separations at which particles can bind to both surfaces. The surface coverage at 
these small separations remains close to 100% for volume fractions (j)* /lO J; <^ 10 0*, while the coverage Coo = 4'/{4' + for 
large separations changes between approximately 9% and 90% in this example. 



fractions ^ 1 of the particles and large binding energies 
U with e'^/"^ ^ 1. Here, d is the particle diameter, and T 
denotes the temperature in energy units. 

Interestingly, the effective adhesion energy (1) of the 
surfaces exhibits a maximum at an optimum bulk volume 
fraction of the particles, see fig. 2. At this volume frac- 
tion, the particle coverage of two planar parallel surfaces 
turns out to be close to 50% for large separations ('half 
coverage'), and 100% ('full coverage') for short, binding 
separations, see fig. 2. Bringing the surfaces from large 
separations within binding separations thus does not 're- 
quire' desorption or adsorption of particles at the optimum 
volume fraction. The existence of an optimum particle 
volume fraction has interesting, experimentally testable 
implications, for example 're-entrant transitions' in which 
surfaces or colloidal objects first bind with increasing con- 
centration of adhesive particles, and unbind again when 
the concentration is further increased beyond the optimum 
concentration. 

Even though our derivations of eq. (1) are based on 
some simplifying assumptions, the effective adhesion en- 
ergy (1) should be applicable in general since it can sim- 
ply be understood as a difference of two Langmuir ad- 
sorption free energies per binding site: (i) the adsorption 
free energy {T /cP) In (l -I- qcjie^^^^) for small surface sep- 
arations at which a particle binds both surfaces with to- 
tal binding energy 2U, and (ii) the adsorption free en- 
ergy {T/cP) In (l -I- qipe^/'^) for large surface separations, 
counted twice in (1) because we have two surfaces. These 
Langmuir adsorption free energies result from a simple 
two-state model in which a particle is either absent (Boltz- 
mann weight 1 — « 1) or present (Boltzmann weights 
q(j)e^^/'^ and q(j)e^^'^, respectively) at a given binding 
site, see e.g. [11]. The factor q depends on the degrees of 
freedom of a single adsorbed particle and has to be deter- 



mined from more detailed binding models. 

In the following, we will study particles of diameter d 
that are attracted to the two surfaces within the binding 
range r. We will first consider a 3-dimensional lattice gas 
model with d = r that leads to relation (1) with g = 1. We 
will then study a more general tube model with arbitrary 
values of d and r from which we obtain another derivation 
of (1) with q = r/d. Finally, we generalize our main result 
(1) to surfaces with specific receptors for solute particles 
or molecules. 

Lattice gas model. — The simplest model for a gas of 
adhesive particles between two parallel and planar surfaces 
is obtained by discretizing the space between the surfaces 
into a cubic lattice. In this model, the hard-core interac- 
tions between the particles are incorporated by choosing a 
lattice spacing d equal to the particle diameter. Each lat- 
tice site then can contain only one particle, and the bulk 
volume fraction of particles is = e^/'^ /{I + e^/'^) where 
/X is the chemical potential. If a particle is located at a 
lattice site adjacent to one of the two surfaces, it gains 
the binding energy U . The separation of the surfaces is 
i — n±d where n± is the number of lattice layers between 
the surfaces. The total number of lattice sites between the 
two surfaces is nynj^. 

The free energy F of this lattice gas of particles between 
the surfaces can be decomposed into a bulk and a surface 
term. The bulk free energy /{,, which is equal to the free 
energy per lattice site in the limit of large nyrij^, has the 
simple form ft = -Tln(l + e^/"^) = Tln(l - 0). The 
surface free energy fg = (F — nj|nj_/f,)/n|| is the excess 
free energy per pair of apposing surface sites in the limit 
of large ri || . The surface free energy depends on the sepa- 
ration of the surfaces, i.e. on the number of lattice layers 
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n±. The surface free energy is given by 
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for nj_ > 2 and by 
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for n_L = 1. The effective adhesion potential of the sur- 
faces, defined as y = fs,n±/(P, thus is constant for n± > 2 
and has an attractive well of depth Ud — {fs.oo — fs,i)/<P 
at nj_ = 1. After replacing the chemical potential /i 
in eqs. (2) and (3) by the bulk volume fraction (j) — 



we obtain 
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For small volume fractions <C 1 and large binding ener- 
gies U with e^/"^ ^ 1, the effective adhesion energy (4) is 
identical with (1) for q ~ I. 

Tube model. — We now consider a semi-continuous 
model to obtain a realistic estimate of the factor q in 
eq. (1) from the separation-dependent effective interaction 
potential of the surfaces. In this model, we discretize the 
space between the surfaces into tubes of the same diame- 
ter d as the particles. We assume that the two apposing 
surfaces are on average parallel and nearly planar, with lo- 
cal curvature radii much larger than the diameter d of the 
adhesive particles. The tubes are oriented perpendicular 
to both surfaces and contain the particles. The length £ of 
the tubes thus corresponds to the local separation of the 
surfaces. The particles can exchange between the tubes 
and with the bulk solution (see fig. 1). We assume that 
the attractive interaction of a particle with the surfaces is 
short-ranged and model this interaction by a square-well 
potential with binding energy U and range r < d/2. A 
particle is thus bound to a surface with binding energy 
U if the separation between the surface and the center 
of the particle is smaller than r + d/2. The pair inter- 
actions between particles are purely repulsive hard-sphere 
interactions. For simplicity, we assume that the tubes are 
arranged on a square lattice. Each tube then occupies a 
volume d'^i. The approximation implied by this discretiza- 
tion should be valid for the experimentally relevant small 
particle bulk volume fractions 4> considered here. Since 
the particle-surface interaction is short-ranged, the gas of 
particles between the surfaces is as dilute as in the bulk 
for large surface separations £ ^ d, except for the sin- 
gle adsorption layers of particles at the surfaces. For large 
binding energies U, these adsorption layers will fully cover 
the surfaces. 

In this model, the problem of determining the effective 
adhesion potential of the surfaces is reduced to calculating 
the partition function of a one-dimensional gas of particles 
in a tube. Since the number of particles in a tube varies. 




Fig. 3: Effective adhesion potential V as a function of the 
surface separation £ for the bulk volume fraction <j!> = 0.01, 
binding energy U = 8T, and binding range r — ci/4 of the 
particles with diameter d. The potential has a minimum at 
the surface separation £ — r + d and attains a constant value 
for separations £ > 2{r + d). The effective adhesion energy of 
the surfaces, Ud, is the difference between the asymptotic and 
minimum value of potential V{£). For large binding energies U 
with e'^^'^ ^ 1, the effective potential has a barrier of height 
Uha at surface separations d + 2r < £ < 2d. 



the suitable statistical ensemble is the grand-canonical en- 
semble in which the chemical potential fi of the particles 
is fixed. The free energy is 
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where -E„ denotes the canonical partition function for a 
system of n particles confined in the tube. The upper 
limit of the sum, \_£/d\ , equals the largest number of hard 
spheres of diameter d that can be placed in the tube of 
length i. The effective adhesion potential of the surfaces 
isV= {F-fbd^£)/d^, where fb = lim^^oo F/{d^£) denotes 
the free energy density in the bulk. The effective adhesion 
potential thus is again defined as the surface contribution 
to the free energy F. With eq. (5), we obtain 
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One-dimensional models for hard spheres have been 
studied extensively [11-16]. Two well-known results that 
we will use in the following are: (i) The hard-sphere gas fu- 
gacity is e'^/'^ « (0 -|- 20^) A/d, up to second order terms 
in the bulk volume fraction [14]. Here, A denotes the 
thermal de Broglie wavelength, (ii) The bulk free energy 
density is ft = —T(f)/{d^{l — (f))), with the denominator de- 
scribing the volume accessible to a sphere [11]. In contrast 
to previous studies, we are interested here in the effective 
adhesion potential V given by eq. (6). With the relations 
(i) and (ii) above, a virial expansion of the potential V 
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Fig. 4: Effective adhesion potential V as & function of the sur- 
face separation I for the particle bulk volume fraction (j) = Q.l, 
binding range r = d/4, and binding energies U = COT (lower 
curve) and U = 0.6 T (upper curve). The global minimum of 
potential V{(.) is located at surface contact £ — for small en- 
ergies U < |r ln(l -I- d/r) (upper curve), and at the separation 
£ = d + r for binding energies U > |rin(l-|-d/r) (lower curve). 



leads to 



V : 



T 



ai 
a2 



(7) 



(8) 



dxi e 



-H{xi)/T 



(9) 



for (f> 1 with the expansion coefficients 
[AZi{£)-£]/d and 

[(2d - £) AZi{e) + K^Z2{1) -id + £^/2] /(f 
The partition function Zi is given by 

with the one-particle configuration energy H{xi) — 
-Ue {d/2 + r-xi) -U<ci{xi-£+ d/2 + r). Here, xi 
denotes the center-of-mass position of the particle in the 
tube, and Q is the Heaviside step function with 0(x) = 1 
for X > Q and 0(a;) — for a; < 0. We have assumed here 
that the tube length £ is larger than the particle diameter 
d. For tube lengths £ > 2d, the partition function Z2 is 

1 j-l-d/2 j-x-2-d 

Z2=^ dx2 da;ie-^(^i'-^)/^ (10) 

^ J3d/2 Jd/2 

with the two-particle configuration energy H(xi,X2) — 
~Ue {d/2 + r- xi) - Ue{x2-e + d/2 + r). For tube 
lengths £ < d and £ < 2d, respectively, the partition func- 
tions Zi and Z2 are zero. 

To perform the integrations in eqs. (9) and (10), it is 
helpful to rewrite the integrand by making use of the sim- 
ple relation e'^®(^)/^ = 1 + (e^/^ - l)e(x). The integra- 
tions then lead to 

p2U/T r 

fn{t-2r) 



(l-e-^/^)f„(^-r) + f„(£) 



(11) 



with the auxiliary function f„(^) = {I — nd)^Q{l — nd)/n\. 
The effective adhesion potential V{t) finally is obtained 
from inserting this result into eqs. (7) and (8). 

For small bulk volume fractions <f> considered here, the 
effective adhesion potential is constant for separations t > 
2(d+ r), with the asymptotic value 
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The adhesion potential has a local minimum 
T [ /r 

d^H' + '^id 

at the surface separation i — r+d, see figs. 3 and 4. At this 
separation, terms of order 0^ are negligible since the two- 
particle partition function Z2 is 0. The virial expansion 
coefficients Oi and 02 then are of the same magnitude, and 
a2(/)^ is much smaller than ai0 for 0^1. 

Depletion interactions are reflected in a second mini- 
mum of the effective potential at surface contact ^ = 0, 
see fig. 4. For surface separations < ^ < d, the potential 
V{£) = £T(j)/[{l — (l))d^] increases linearly with £ due to 
the depletion forces. By definition, the effective potential 
vanishes at surface contact ^ = 0. The two minima thus 
have equal depths for V,nin = 0, i.e. at the binding en- 
ergy Uq « ^Tln(l -I- d/r). For binding energies U > Uq, 
the global minimum of the potential V{£) is located at 
i = r + d. The effective adhesion energy of the surfaces 
then is Ud — Voc — Knin with Voc and Vniin as in eqs. (12) 
and (13). For large binding energies U with e'^^'^ ^ 1, we 
obtain the effective adhesion energy (1) with q = r/d. 

The adhesion energy Ud is maximal at the bulk volume 
fraction (j) — (j)* with (j)* « e~^/^ d/r, see fig. 2. At this 
optimal volume fraction, the coverage of the unbound sur- 
faces. Coo = -{d'^/2){dV^/dU) « (j)/{cj) + (j)*), is 50%. In 
contrast, the coverage of the bound surfaces at separation 
£ = d+r is c„,i„ = -(dV2)(9Knin/a(7) « cj)/{^+c^*e-^'^) 
and approaches 100% for (j) fa (fi*. 

Besides the two minima at ^ = and £ = d+r, the effec- 
tive adhesion potential has a barrier at surface separations 
d+2r < £ < 2{d+r) for large binding energies, see fig. 3. 
At these separations, only a single particle fits between 
the surfaces, but this particle can just bind one of the sur- 
faces. The particle thus 'blocks' the binding site at the 
apposing surface, see particles in the center of fig. 1. At 
the surface separation £ = 2d within the potential barrier, 
the effective adhesion potential V{£) attains the value 

lH'^K^-d{' 

for small bulk volume fractions <j). At large binding ener- 
gies U with e^/^ ^ 1, the barrier height C/ba = Vco — Vha 
then is 

r {l + ^e^/^r/df 
d2 ^ l + 2<j)eU/Tr/d 



(14) 



(15) 
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Fig. 5: Two surfaces, e.g. lipid membranes, with specific re- 
ceptors for solute particles or molecules in the bulk. For small 
surface separations, a particle can link the surfaces together by 
binding to two receptors at apposing surface sites (particle on 
the left). 



Surfaces with specific receptors. The effective 
adhesion energy (1) can be generahzed to cases in which 
the adhesive particles specifically bind to receptor sites or 
molecules at the surfaces, e.g. to receptors in fluid lipid 
membranes. An example are biotinylated lipids, which 
can be crosslinked by streptavidin [f7, f8]. In principle, 
the two membranes can contain the same type of recep- 
tors, as in fig. 5, or different types of receptors. We charac- 
terize the receptors in membrane 1 by their binding energy 
Ui and chemical potential vi, and the receptors in mem- 
brane 2 by the binding energy U2 and chemical potential 
1^2- For simplicity, we assume the same binding range r 
for both types of receptors. The chemical potential of 
the receptors is the free energy difference between a mem- 
brane patch of size d? containing a receptor molecule, and 
a membrane patch of the same size without receptor [f9]. 
fn the absence of adhesive particles, the chemical poten- 
tials i^i and 1^2 are directly related to the area fractions 
ai = e''i/'^/(f + e"^/^) and = e''^/'^/(f + e''^/^) of the 
receptors in the two membranes. 

Let us first consider the Langmuir adsorption free 
energy of a single membrane with receptors. Since 
each patch of the membrane can attain four possible 
states, the Langmuir free energy per patch area (P is 
(T/d2)ln(l + e'^'/^ + g0 + q0e('^'+''')/^) with i = 1 or 
2. Here, e''*/"^ is the Boltzmann weight for the state in 
which a receptor is present in the patch but no particle 
is within binding range, q(j) is the Boltzmann weight for 
having a particle within binding range in absence of a re- 
ceptor, and 6^'^'+'^*-'/"^ is the Boltzmann weight for a 
receptor bound to an adhesive particle. The prefactor q 
depends on the degrees of freedom of bound particles. In 
our model, we obtain q = r/d, see previous section. 

At large separations, the Langmuir adsorption free en- 
ergy of two membranes simply is the sum of the adsorp- 
tion energies above. At the optimum binding separation, 
however, the Langmuir adsorption energy of two apposing 
membrane patches is (T/d2)ln[(l + e''i/^)(f -I- e"^/'^) + 
g0(l + e('^i+'^i)/^)(f + e('^2-(-i'2)/T)]^ before, the effec- 
tive particle-mediated adhesion energy of the membranes 
is the difference between the Langmuir adsorption energies 



at large separation and the optimum binding separation. 
In the symmetric case with binding energies U — Ui = U2 
and chemical potentials v = vi = U2, '^e obtain 

T l + <i>{l-a + ae^'Tf r/d 

Ucf w — In ^ 7j (16) 

{l + ^{l-a + aeU/T)r/d) 

with a = e'^/'^ /{I + e'^/'^). The effective adhesion energy 
(16) can be obtained directly from eq. (1) by replacing the 
Boltzmann factor e^^'^ with the term 1 — a -I- ae^/'^ . For 
a = \ where each membrane patch contains a receptor, the 
eqs. (1) and (16) become identical. The effective adhesion 
energy (16) is maximal at the bulk volume fraction (p* w 
(d/r)/(l-a + ae^/^). 

Discussion and outlook. — We have determined the 
effective interactions between two surfaces in contact with 
adhesive particles. Our results apply to a wide range of 
surfaces, which may be soft or rigid, and planar or non- 
planar. For non-planar surfaces, the local surface sep- 
aration t varies, and the total particle-mediated interac- 
tion energy is obtained by integration over the separation- 
dependent local interaction V{t) defined in (6). The over- 
all interaction potential of the surfaces is a superposition 
of the particle-mediated interactions with direct interac- 
tions, such as van der Waals, electrostatic and hydration 
forces. For large particle radii d, the particle-mediated in- 
teractions should dominate since the particle-bound sur- 
faces have a separation close to d-l-r, where r is the range 
of the attractive surface-particle interaction. The interac- 
tions induced by the particles then can be measured di- 
rectly, e.g. via the surface- force apparatus [7], or inferred 
from the phase behavior of colloidal systems [5,6]. 

The effective adhesion energy (1) has been obtained 
from two different models: (i) from a 3-dimensional lat- 
tice gas model; and (ii) from a more elaborate tube model 
for hard spheres. In both cases, we discretized the space 
in order to incorporate the hard-core interaction between 
the particles in an analytically tractable manner. In the 
lattice gas, we discretized all three coordinates of the par- 
ticles. In the tube model, the coordinate perpendicular 
to the two surfaces was taken to be continuous. As far 
as relation (1) is concerned, the only difference between 
the lattice gas model and the tube model is that they give 
somewhat different expressions for the dimensionless co- 
efficient q. Thus, it is rather plausible that the effective 
attractive interaction (1) also applies to hard spheres with 
three continuous spatial coordinates. This proposition can 
be checked by Monte Carlo or Molecular Dynamics sim- 
ulations, which have been previously used to study hard 
spheres confined between nonadhesive surfaces [20-24]. 

We have neglected that flexible surfaces such as lipid 
membranes can wrap around adhesive particles [25-28]. 
A partial wrapping leads to effective, surface-mediated in- 
teractions between the adsorbed particles [29-31], and to 
cooperativity effects in adsorption [32]. We have also as- 
sumed that the receptors considered in the last section 



P-5 



Bartosz Rozycki, Reinhard Lipowsky and Thomas R. Weikl 



are rigid molecules and not flexible tethers [33-35], which 
seems to be a good assumption for most biological recep- 
tors. In addition, wc have neglected direct, long-ranged 
interactions of the adhesive particles, e.g. electrostatic re- 
pulsion of charged particles. At small bulk concentrations 
of the particles, repulsive interactions will mainly affect 
the packing density of the adsorbed particles and, thus, 
the concentration of available binding sites. For charged 
solutes, such as multivalent ions [8] or charged proteins [36] 
adhering to lipid membranes, the average particle separa- 
tion at maximum surface coverage is affected by the salt- 
dependent screening length of the solutes at the surfaces. 

Thermal shape fluctuations of lipid membranes lead to 
an unbinding transition if the fluctuation-induced repul- 
sion exceeds the effective adhesion energy Ugf. The charac- 
ter of the transition depends on the barrier in the effective 
adhesion potential. According to scaling arguments [37], 
the unbinding transition of the membranes is discontinu- 
ous for strong barriers with J/ba^ba > cT^/k where c is a 
dimensionless coefficient of order 0.01 [38], and continuous 
for weak potential barriers. The thickness of the barrier 
here is ?ba ~ d, see fig. 3, and k is the bending rigidity 
of the lipid membranes. At the optimum bulk volume 
fraction (f)* e^^^/"^ d/r, the effective barrier strength is 
C^ba<^^ ~ 0.3T and thus clearly beyond the threshold for 
discontinuous unbinding, since the bending rigidity of lipid 
membranes typically is between 10 and 20 T. 

In this letter, we have considered equilibrium aspects 
of adhesion. The unbinding dynamics of surfaces with 
multiple receptor-ligand bonds under a pulling forc;e has 
been studied in [39-41] . If receptors in apposing surfaces 
are connected via adhesive solutes, each bond between the 
surfaces consists of two molec;ular bonds in series. Such 
serial bonds can break at either of their bonds, and there- 
fore have been found to break earlier than single receptor- 
ligand bonds under an applied force [9,10]. 
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